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Abstract. We consider discrete Schrodinger operator J with 
Wigner-von Neumann potential not belonging to . We find asymp- 
totics of orthonormal polynomials associated to J . We prove the 
Weyl-Titchmarsh type formula, which relates the spectral density 
of to a coefficient in asymptotics of orthonormal polynomials. 



1. Introduction 

In recent papers on Jacobi matrices |I] , [E] , [13] , [25] , [2l] , [22] new re- 
sults were found on the asymptotics of generalized eigenvectors of these 
operators. For real sequences {anjj^i and {6„}^^ the Jacobi operator 
J = J(a„, 6„) is defined in the Hilbert space by the formula 

As it is well known the spectral analysis of Jacobi operators is strongly 
related to the study of the asymptotics of generalized eigenvectors. In 
this work we concentrate on the discrete Schrodinger operator J = 
J^{1, bn) with the Wigner-von Neumann potential 



, , , csin(2a;n + 6) 



(2) 



7 e Q; , 2u; ^ ttZ and {gj^^^ G l\ 



where c, u, S, g„ G M. JT" is the Jacobi operator given by 

/gx {Ju)i = biUi + U2, 

{Ju)n = Un-l + bnUn + M„+i, 71 > 2. 



1991 Mathematics Subject Classification. 47B36,34E10. 

Key words and phrases. Jacobi matrices, Asymptotics of generalized eigenvec- 
tors, Orthogonal polynomials, Weyl-Titchmarsh theory. Discrete Schrodinger oper- 
ator, Wigner-von Neumann potential. 



2 



JAN JANAS AND SERGEY SIMONOV 



and has a matrix representation in the canonical basis {e„} 

\ 

1 h 1 ■■■ 

J 



oo 
n=l 



of f 



1 





1 

1 





1 

&3 



V 



Since J7 is a compact perturbation of the free discrete Schrodinger op- 
erator its essential spectrum is the interval [—2; 2]. Frequency 2u in the 
potential produces (in general) four critical (or resonance) points inside 
this interval: ±2 cosw, ±2 cos 20;. At these points the resonance occurs 
and the asymptotics of the generalized eigenvectors changes (analogous 
phenomenon in the continuous case is very well studied, see [lO], [8], 
Chapter 4 and [21], Theorem 5) and an eigenvalue can appear under 
certain additional conditions. 

In the present paper we are interested in asymptotics of orthogonal 
polynomials P„(A) associated to J', which are 

Pi(A) := l,P2(A):=A-6i, 

P„+l(A) := (A - bn)Pn{X) - Pn-l{X), ^ > 2, 

and the relation of this asymptotics to the spectral density p'(A) of J'. 

Our main result (see Theorem [1] on page [H] and Theorem [2] on page 
20] for the exact formulation) is that there exists a function F ( the 
Jost function ) such that for a. a. A G (—2; 2), 



(4) 



and 



P'(A) 



2vr|F(z)| 



where 



zF{z) 1 



z"^ z' 



exp 



1-27 



27 



+ ^^2;" exp 
z^ — \ 



li2\z)n 



1-27 



A = 2 + 



A 



1-27 



+ 0(1) as — > 00, 



and 

(5) 



c^z^{l + z^)e 



2iuj 



4{1 - z^){z^ - e^'^){l 



This result is an analog of the classical Weyl-Titchmarsh (or Kodaira) 
formula for the differential Schrodinger operator on the half-line with 
summable potential (see [28], Chapter 5 and [T7]). 
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We consider 7 G (|; |) for the following reasons. If 7 > |, then 
{bn}^=i G and the sum 

00 

n=l 

converges. This situation was studied using a completely different 
method by Damanik and Simon in [7] for a class of more general Ja- 
cobi matrices, with arbitrary sequence {bn}'^=i G for which the se- 
ries Yl'^=i is convergent. In this case the asymptotics of polynomials 
P„(A) is simpler (of Szego type): 

„ /^^ zFiz) 1 zFiz) „ , , 
(6) P„(A) = ^ ■ - + ■ z- + oil) as n ^ 00 

, and this was proven in [7] not for a. a. A G (—2; 2), but in the sense 
of the convergence in L2((— 2; 2), p'{x)dx) (Theorem 8.1). If one proves 
that (E]) holds for a. a. A G (—2; 2) (which for the potential of the form 
(II]) is much simpler than the analysis that we develop in the present 
paper), then (jl]) can be deduced from results of Damanik and Simon 
(Theorem 5.6 and Theorem 8.1) in a non-trivial way (this was pointed 
out by Dr. Roman Romanov in private communication). On the other 
hand, the type of asymptotics for 7 G (|; 1] is the same as in the simple 
case 7 > 1. For 7 = | the type of asymptotics changes. This happens 
also for every 7 = j, I G N\{1}. We are forced to consider different 
(depending on /) number of terms in the asymptotic expansion that 
we use, see ( fT5l) . The greater / is, the more terms are significant. The 
method that we use works for every 7 > 0, but we restrict ourselves to 
the case 7 G (|; |) to show how it works in the general case. However, 
in the final section we state the corresponding result for 7 G (|; 1] 
(Theorem [3]) and indicate how to simplify the proof for that case. The 
main idea of the method is inspired by |6] and uses the discrete version 
of a change of variables introduced by Eastham in |8]. 

We could as well consider a finite sum of terms like ([T]) as the poten- 
tial, 

, Q sm{2uin + 61) , 

bn = }_^ + 

1=1 

with the same conditions imposed on ci,ijJi,'yi,6i and This 
would increase the number of critical points and complicate the nota- 
tion and calculations, so we restrict ourselves to the case of one such 
term only. 

In the continuous case (differential Schrodinger operator on the half- 
line) Wigner-von Neumann potentials were studied in numerous works: 



4 



JAN JANAS AND SERGEY SIMONOV 



, [2], [3], [n], [IE], [19], 120], [6], [21], and formulas for the spectral den- 
sity analogous to (jl]) were obtained in different variations in [22], [2] and 

i- 

The second name author plans to use the Weyl-Titchmarsh type 
formula to study the spectral density of the Jacobi matrix J . For the 
case 7 = 1 there are only two critical points: ±2 cos a;. Analyzing the 
change of asymptotics of Pn(A) as A approaches the critical value Aq G 
{±2cosa;} (remind that the type of asymptotics changes at A = Aq), 
we plan to prove in the forthcoming paper for the case 7 = 1 that 

|e| 

(7) p'l-^) ~ const • I A — Aol^i"'""^! as A — t- Aq, 

if {P„(Ao)}^i is not a subordinate solution of the spectral equation. 

Genralization to the case 7 < 1 is possible, but the problem seems 
to be much more difficult than for 7 = 1. 

2. Preliminaries 
For every complex A the spectral equation for J 

(8) Un-\ + hnUn + Mn+1 = Am„, n > 2 

has solutions Pn(A) (orthogonal polynomials of the first kind) and 
Qn(A) (orthogonal polynomials of the second kind) such that 

Pi(A) = l, P2(A) = A-6i 
Qi(A) = 0, g2(A) = i. 

For non-real values of A there exists m(A) (the Weyl function) such 
that 

g„(A) + m(A)P„(A) 

is [T] the unique (up to multiplication by a constant) solution of dH]) 
that belongs to The Weyl function and the spectral density of J 
are related by the following equalities: 

^ r dpix^^ ^ ^ 

J X- X 

p'(A) = -Imm(A + iO) for a.a. A G M. 

71 

For every two solutions u and f of ([8]) with the same parameter A the 
(discrete) Wronskian defined by 

W{U, V) := UnVn+l - Un+lVn 

is independent of n. For polynomials P(A) and Q{\) the Wronskian is 
equal to one for every A. 



weyl-titchmarsh type formula 
3. Reduction of the spectral equation to a 

SYSTEM OF the L-DIAGONAL FORM 

Consider the spectral equation for J , 
(9) + 6„,M„ + = Am„, n > 2. 

Let us write it in the vector form, 



(10) 



1 

-1 \-hr. 



Un-1 
Ur, 



n>2. 



Consider a new parameter 2; G D (we denote by D the open unit disc) 
such that 

\ = z + - 

z 

and conversely 



where the branch of the square root is chosen so that 2; G © for A G 
C \[-2; 2], i.e., ^ = -i for A = 0. Let 



z'-l \ - 



which is equivalent to 



Un+1 



1 1 

1 ~ 



Un 
Un+1 



Vn, n>l. 



We make this substitution to diagonalize the constant part of the co- 
efficient matrix 

-1 



1 
-1 A 

Equation (ITU]) becomes 

i 

z 

z 



1 1 

1 ~ 



11 



i 

z 

z 



1 z2 
-1 -z' 



1 1 

i z 



Vn.n > 1. 



The goal of the present section is to reduce the system flTT]) to the 
" L-diagonal form" . If we put 

Wn := T~'^{z)Vn, 



then (jTTj) becomes 

Wn+l = T-^liz) 



i 

z 



+ 



Z^ - 1 



1 z' 



Tn{z)Wn- 
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The system is in L-diagonal form if the coefficient matrix is a sum of 
diagonal and summable matrices. So we have to find matrices Tn{z) to 
provide this property. This is possible not for every 2; 7^ 0, 1, —1. 
Let us denote 

(12) U :=€ \{0, 1, -1, e^''^, -e^'^, e^^'"", -e^^i^}. 

Lemma 1. Let z G U. For every n G N there exist matrices Rn\z) 
and invertible matrices T„(z) such that 

T^{z),T~\z) = I + o{l), 



(13) 

as n — 7- 00 and 



i 

z 



+ 



Z^ - 1 



1 I f^2{z) 

J- -r „27 




Tn{z) 





z I 



M2(^) 



,27 



where jj,2{z) is given by ([5]). Tn{z) and Rn\z) are also analytic in U for 
all n, and on every compact subset of U the estimate f lT3|) is uniform 
with respect to z. 

Proof. Let us denote 

A{z) 

N2{z) := 

N^2{z) :-- 
?(o) 



i 

z 

z 

^gi{2aj+5) 
2i(22_l) 

gg-i{2aj + 5) 

^ ~ 2j(z2-l) 

1 



1 ^2 
-1 

1 

-1 



^^-1 \ -1 -z' 



so that 



1 



-2iLon 



N2{z) + ^--N_2{z) + R'^\z). 



We will find T„(z) in two steps. 
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At the first (and main) step of the construction of Tn{z) let us find 
matrices Tn\z) such that 



(14) [Ti'U^ 



i 

2 

z 



+ 



bn+l 

- 1 



1 z' 
-1 -z^ 



Ti'\z) 

n M + ^+i?(^)(z) 



where V{z) does not depend on n and Rn\z) is summable. 

Following the ideas of [S] and [0] we look for Ti^'' (2;) of the form 



(15) r«(^):=exp 



p2iu)n p—2ium 

X2{Z) + X_2(^) 



Hum 



+ ^X4(Z) + 



-X_4(Z^ 



where X±2{z) and X-|-4(z) are to be determined. Define 
1 



(16) M±4 := -(AX±2' + e±^^-X±2'A) + N±2X±2 



±2iw ■ 



and 



(17) V := -(A(X2X_2 + X_2X2) + (X2X_2 + X_2X2)A) + X2X„2 

+X_2X2-(e2*"X2X_2+e-''"X_2X2)-(e2*"X2AX_2+e-''"X_2AX2). 

Take expansions of Tn^^ and (T^^\)~^ as n — )■ 00 up to the terms of the 
order After a long but transparent calculation we have: 



(18) (ri^i 



-1 



^2iu)n 



-2iu)n 



A + X2 + 



X_2 + 



^,2iujn 



~2 
-2ium 



= A + [Xa + AX2 - e^*^X2Al ^ 

4ja;n — 4jajn 

+ ^[M4 + AX4 - e^-X4A] + -- 

^/7 ^^7 

+ -^ + 



-[X_2 + AX_2-e-'*"X_2A] 
-[M„4 + AX_4 - e-^^'^X_4A] 













as n — 00 



since = 0(|g„+i|). 

We want to cancel the coefficients at e^^*"^" and e^"^*"^" in (fTSll by 
suitable choice of X-|-2(-2) and X-|-4(z), respectively. To this end, four 
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conditions should be satisfied: 



gSi^XaA - AX2 



No 



-2iui 

e^'^X^K - AX4 = 
-^'"^X^k - AX 



M4, 

e"''"^A_4A - AA_4 = M_4, 

We use the following lemma to solve them. 
Lemma 2. If fi ^ 1, z"^, ^, then the matrix 



X 




satisfies the equation 



/xXA - AX 



2^^— 1 
.122 

fll fl2 
/2I /22 



Proof. The assertion can be verified by direct substitution. 
It follows that we can take 



□ 



(19) 



and 



X2iz) 
X.2iz) 



„i{2ui+6) 



'1 



J2iui_^2 
-«(2"+«) / T=2Er 



-2iu 



-2iijj _ ^ 



X. 



z{M±4{z))ii z{M±4,(z))i2 

iM±4{z))22 

z2 2(e±4i^_l) 



z{M±4iz))21 



where (M^ (2;)) 1112.21,22 are the entries of the matrix M±i{z), which 
are given by (ITBi) and (IT^ . As we see, X-|-2(-2) are defined and analytic 
in C \{1, —1, e^*"^, — e^*"^}, M±i{z) and V{z) are defined and analytic in 
C \{0, 1, -1, e^'^, -e^'^} and X±4(^), i?f (z), Ti'^(z), {Ti,^\z))-^ are de- 
fined and analytic in U . 

With this choice of Tn^ iz) the remainder 



J2iu}n 



K{z) + — -X2(^) 



-2iijjn 



+ 



-X_2(z)+4°)(^) 



T«(.)-A(.)- 



27 



satisfies the estimate 

= O (^-^ + |g„+i|^ as n ^ 00 
uniformly with respect to z on every compact subset of U . 
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At the second step of the construction of matrices Tn{z) let us con- 
sider 



z) := exp 



n 



27 



with some Y{z), which is to be determined. Taking the expansions of 

(2) ^ 



Tn^^ and (T!^+i) up to the order of -57, we have: 



n 



27 









- Wn+llj 



Let us cancel the anti-diagonal entries of — [Y, A] by the choice of Y. 
This leads to the equation 



We can take 

Y(z) 



Vu 
V21 



V,2(z 



z^-l\ -V2i{z) 

What rests is 

/ 1 

diag V{z) = fi2{z) ( ^ ^ 

which can be seen from f ll7l) and fll9p by a straightforward calculation. 

Matrices Y{z),Tn^\z) and {Ti^\z))~^ are defined and analytic in U. 
The remainder in the system after the transformation, 

R^^\^) ■■= (a(.) + ^ + T(^)(.) 

diag\^(2;) 



-A(z)- 



^27 ' 
satisfies the estimate 

\\R^n\z)\\ = O (^-^ + |g„+i|^ as n ^ 00 

uniformly with respect to z on every compact subset of U. Taking 

Uz) := n'\z)T('\z) 
we complete the proof. □ 
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Finally, we have come to the following system of the L-diagonal form: 



(20) 



1 + 



,27 











z 1- 



It is easy to check that for z ETdU the value fi2{z) is pure imaginary. 

4. Asymptotic results 

In this section we prove several results needed for the analysis of 
the system fl20l) . They are more or less standard, and the approach 
is similar to [12], [26] and [1]. In the cited papers, the existence of a 
base of solutions with special asymptotic behavior is proven. Here we 
find the asymptotic of (roughly speaking) generic solution defined by 
its initial value. 

Let us use the following notation. 



E :=o, n 

n=l n=l 

and for rii, n2 > 1, 

712 n2 ni — 1 

n=ni n=l n=l 



n2 

n 

n=ni 



n2 



n ( n 

n=l \ n=l 



-1 



The first lemma is a kind of discrete variation of parameters. 

Lemma 3. Let / G and let the matrices A„ be invertihle for every 
n>l. If for every n > 1 



(21) 



n— 1 / n— 1 



k=l \l=k+l 



then 
(22) 

for every n>l. 
Proof. Consider 



{An + R 



,(1) 



We can rewrite fl22l) as 
(23) vUl 



ffl-l 

Ha, 



Ha. 



-1 



Xr. 



Rn-^n- 
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At the same time ( 12T1) is equivalent to 



11 



n—l / k 



(24) 



k=l \l=l 



Clearly, ([23]) follows from (EH). 



□ 



Remark 1. Lemma\^ says that x given by fl2Tl) is the solution of the 
system fl22|) wi/i xi = f . Every solution x of (l22l) can 6e represented 
in the form fl2T|) wt/i f := Xi. 



In what follows let us consider systems with matrices A„ of the form 

Ar, := 



A„ 



where A„ G C. The following lemma gives an estimate on growth of 
solutions of the system (|22|) . 



Lemma 4. Lei 



°° II R II 



fc=l 



and let there exist M such that for every m > n 

m 



l=n+l 

Every solution x of the system 

A„ 

An 

5' estimate: 

n 



(25) Xn+l = ( 'n" 1 ) + -Rn 

satisfies the following estimate 

(26) ||x„|| < (H exp I (1 + M^) ^ 
Proof. Let 



Xn, n > 1 



fc=i 



:i + M2)||xi||. 



/ 



/i 

/2 



Define 



n-1 

n A. 
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Then {?/„ } satisfies the equation (by using Lemma [3]) 

/I \ n-l / 1 \ 

Consider this as an equation in the Banach space /°°(C^). Denote 

/i 



the vector 



n-l 



n=l 

1 



and the operator V : {u„}~=i <( E I g J- I fr^^ 



3? 



l=k+l '■ 

Equation ( 127|) reads in this notation: 

The powers of the operator V can be estimated as follows: 



n=l 



oo 



ml 

and so (/ — V)^^ exists and 



m 



\\{I-V)-'\\,^<exp[{l + AP)J2 

Therefore 



oo ,, „ 



A;=l 



lA 



< exp (^(1 + M') ^ j (1 + 

Returning to the solution x we have arrived at the desired estimate 

□ 



The following lemma gives asymtotics of the solutions of the system 
Lemma 5. Let 

k=l ' 

and let there exist M such that for every m > n, 

m _ 

(28) n i^'i £ M- 



l=n+l 
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Suppose that x is a solution of the system 



(29) 

a) If 

(30) 

then 



K 



+ Rn 



Xn, n > 1. 



< oo, 



-1 



-1 



lim TTA 



l\ Xn = Xi 



RkXk 



k=i \i=i 

(the limit and the sum both exist and the equality holds). 
b)If 



(31) 
then 



JJ Ai = oo, 



1 . 

lim — - — 

n 

1=1 



1 




Xi 



fc=i 



RkXk 

~k 



(^t/ie /zmzi and the sum both exist and the equality holds). 
Proof. Case a). Equation f l24|) can be rewritten as follows: 

(n-l \ "1 n-l / fc \ ^1 

Y[Ai] Xn = Xi + Y,\Yl^l] RkXk. 
1=1 / k=l \l=l / 

Let us show that the sum on the right-hand side is convergent. By 
Lemma HI 



Ha, 

J=l 



RkXk 



< 



\\Rk 



n 

j=i 



A; 



X 



exp (l + M^)^^ 



11^. 



1 |A. 



(1 + M^)||a:i||. 



Since 



n 

J=l 



1 

flAf 

1=1 



< 



\ 



j=i 
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is bounded by hypothesis, we have: 



Ha. 



RkXk 



< const 



\Rk 



\Xk 



which is summable. Therefore the hmit in (1321) as n — )■ oo exists. 

Case b). Consider the sum on the right-hand side of ( l27j) . Lemma H] 
yields: 



1 




n— 1 

n ^ i 

=fc+i ' / 


Rk 


< 


\\Rk 


\Xk\ 



(2) 

Vk 



< (1 + 



\Rk\ 



I (2) I 

IVk \ 



\Rr. 



m=l 



\X„ 



\Xi I 



which is summable. Since fl3T|) holds, by the Lebesgue's dominated 
convergence theorem there exists the limit as n — )• oo in fl2l 

oo 



lim yl; 



(2) 



1 




Xi 



+5: 

k=l 



1 




Rk (2) 

-^Vk ■ 



Returning to x from y*^^ we obtain the assertion of the lemma in the 
case b). □ 

Remark 2. Condition ( 128|) together with (130|) or (13T]) is a case of the 
standard dichotomy (Levinson) condition, cf. [4j, [12j, [26], [27]. 

5. ASYMPTOTICS OF POLYNOMIALS, JOST FUNCTION AND THE 

SPECTRAL DENSITY 

In this section we apply the results of the previous two sections to 
find asymptotics of polynomials Pn{X) associated to the matrix J' and 
to prove the Weyl-Titchmarsh type formula for the spectral density. 

Theorem 1. Let be given by and the condition ([2]) hold. Let 
Pn{X) be orthonormal polynomials associated to the Jacobi matrix J 
given by ([3]) . Then for every z E 3 H U ( where U is given by (|T2|) ) 
there exists F{z) (the J ost function) such that: 

• if ze D\{0} (i.e., X = z + ^ e C\[-2;2];, then 



(33) Pn 



1 

z + - 

z 



zF{z) 1 



exp 



1-27 



+ 



exp 



1 - 27 
Re^2{z)n^''^^ 
1 - 27 



as n —7- 00 
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(where fi2{z) is given by 1^), 

if z eJnU (i.e., \ = z+\e (-2;2)\{±2cosw,±2cos2a;};, 
then 



(34) P„(. + i)=i^lexp"'='^'« 

' z I I — Z^ Z"' 



1-27 



+ 



zF{z) 



1 - 27 



oil) as — )■ 00. 



Function F is analytic in and continuous inT r\U . 

Proof. For z & U every solution u of the spectral equation (Q corre- 
sponds to the solution w of f l20|) by the equality 



Un 
Un+l 



1 ^ ) Tn{z)Wn. 



Let us define the solution ^9(2;) of ( l20l) that corresponds to polynomials 
Pn(A): 



(35) 
Define 



(36) Xn{z) := 
and 

R^^\z) ■.= R^^\z) + 
so that the system (l20l) reads: 



1 1 v7 ^4^+^) 

^n+l + l) 



i exp (f$J((n + 1)1-27 _ ^1-27)^ ^ if n > 2, 
^exp(f^2i-27V if ^ = 1 



z 1 







(37) 



Wn+l 



Xn{z) 







An(2) 



^27 



w„, n > 2. 



Let us check that Lemmas H] and [5] are applicable to this system. With 
the definition fl5B]) . the product of diagonal entries looks simple: 



n-l 



1 ffi2{z)n'-^^ 
exp 



1=1 



z" 



1 - 27 



n > 2. 
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Let K he a compact subset of D fl t/. We have for m >n: 



where we used the inequahty 

(m + 1)1-2-^ - (n + 1)1-2^ < (m - n)^-27 

(which holds because the function x i— )> x^~'^'^ is concave). Further, 
since Reii2{z) = for 2; e T, the function 

2; H- 



1 - 1^1 

is smooth on U, hence there exists Ci{K) such that for every z & K 

\RepL2{z)\<c,{K){l-\z\). 

Also for every z 

\z\ < el"l-\ 

therefore 

n |A,(^)|>exp (1 - \z\) (m - n - ^^{m - n)'-'-^^ 



l=n+l 

Let 



C2iK) :=sup (^x'-^^-x), 



which is finite. We have: for every z & K and m >n, 

(1 - 1^1) (m-n- - ^)'"''') ^ H2(^); +00) 

and 

m 
l=n+l 

Further, 

\\R^^\z)\\ - O + as n ^ 00 
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1 (^1 + /laM 

Z V 



A„(2;) 



1 H — exp 



27 



1 



^27+1 



as n — > oo 



Analogously 



and finally 



^ _ ^2(^) 



o 



27+1 / ' 



ll^fWII=o(^ + |g„+i|) 



n 



as ri — )■ oo 



uniformly with respect to z & K. Also for every z E K and n, 



|A„(^)| 



< |z| exp 



2|Re/X2(^) 
1 - 27 



Thus the sum 



E 



n=l 



as a function of z is bounded on K (in fact it is continuous in U). For 

zeTnK 



1=1 



is bounded, while for z ^3(1 K 



— > 00 as n — )■ cxD. 



1=1 

We see that Lemma |5] is applicable. It yields: for every z E K there 
exists 



(38) $(.):= [ ^,{z) + ^."exp (- ^Y^;^ j R^^\z)^^{z 



n=l 



and for every z eT (1 K there exists 



$(2;) := I ^1(2;) + ^ — exp 



IJ,2[Z)n 



1-27 



1 - 27 
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and solution (p{z) has the following asymptotics. For z E T n K (case 

(a) of Lemma E]), 

(39) 



,nexp(-H^(g^i^) j vv-^i^; 

as n — )■ oo and for z eBXI K (case (b) of Lemma E]), 

(40) M.) ^ lexp (^^4^>^) (( ) + 0(1)) a, „ ^ oo. 

As we see, Lemma H] is also applicable. Let 



C4(K) :=exp((l + e2^^(^)) C3(ir)) (l + e^'^^^^^) max ||</'i(;2)||. 



Lemma H] yields: for every z E K and n 

II f \\\ ^ ^4(^) /Re/i2(z)n^-27 
||¥'n(2:)|| < ^— exp 



1 - 27 

Consider the expression for ^{z), fl38|) . We have: 

,1-27 

z exp 



27 



<C4(ir)||i?(f)(^)ll 



O { + \qn+i\ ] as n oo 



uniformly with respect to z E K. It follows that the function $ is 
analytic in the interior of K and continuous in K. Since the set K C 
D n f/ is arbitrary, $ exists and is continuous in © fl [/ and is analytic 
in ©\{0}. Asymptotics fl39l) holds for every z G Tflf/ and asymptotics 
(140|) holds for every ^ G ©\{0}. 
If follows that for z G ©UO}, 



^exp — 1 ^ ]Tn{z 



z 

z 



as n — )■ oo, since T„(2;) = / + o(l). If we define 

(41) F{z) := <^(^)^^, 

then we arrive at the first assertion of the theorem. 
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It also follows in an analogous fashion that for z & T (lU, 
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(42) 



{Z + i) 



+ $(2)2" exp 



1-27 



^■2{z)n 



1-27 



1-27 

1 

z 



+ 0(1) as n — )• 00. 



1-27 

The first component of this vector equality describes the asymptotic 
of Pn {z + j), and to complete the proof we need only the following 
lemma. 

Lemma 6. For every z eT r\U , 



$(^) = ^[z). 

Proof. This follows from the fact that values of polynomials Pn {z + ^) 
for 2; G T are real. Consider the imaginary part of the first component 
of 



^{z)-^{z) 1 
= exp 

2i z-^ 



fi2[z)n 



1-27 



1 - 27 



+ \. — —z"" exp 



2i 



fJ'2[z)n 



1-27 



1 - 27 



+ 0(1) as 77, — )■ cxD. 



Suppose that 



Then 



2;^"' exp ( 2iarg ( - 



$(2) ^ $(^). 



— )■ 1 as n — )• 00. 



Let 



We have: 
As well, 



$(z) := exp 2iarg (^$(2;) - $(2' 

Kz) := 

^ 1-27 

|,_2;2ngA"-^"^^ 1 as 77, 00. 



1 

^27 



It follows that z'^ = 1, which is a contradiction. Therefore 

^(z) = Wz). 



□ 
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This completes the proof of the theorem. 



□ 



The following (final) theorem gives a formula of the Weyl-Titchmarsh 
(or Kodaira) type for the spectral density. It follows from asymptotics 
of orthogonal polynomials given by Theorem [1] in a standard way (see 
[28] , Chapter 5, and [17J) and contains the Jost function, which appears 
in the expression for the asymptotics of Pn(A). 

Theorem 2. Let be given by ([1]) and the condition (|2]) hold. Then 
the spectrum of the Jacobi matrix J given by ([3]) is purely absolutely 
continuous on (— 2; 2)\{±2 cosw, ±2 cos 2ci;}, and for a. a. A G (—2; 2) 
the spectral density of J' equals: 



(43) 



P'(A) 



2ti 



2 ' 2 



(Weyl-Titchmarsh type formula), where the Jost function F is de- 
fined in Theorem [H The denominator in (H3|l does not vanish for 
A G (-2;2)\{±2cosu;,±2cos2w}. 

Proof. If we rewrite ( l43l) in terms of the variable z, it reads: 

1-z^ 



(44) 



P 



1 

z + - 

z 



2TTiz\F{z)\'^ 



Polynomials of the second kind have asymptotics of the same type as 
polynomials of the first kind. Indeed, define the cropped Jacobi matrix 
J'l as the original matrix J' with the first row and the first column 
removed, 

/ 62 1 ■ ■ ■ \ 
1 63 1 ■■■ 

1 64 ■■■ 



Jl 



\ ■ ■ ■ ■■■ J 

Polynomials of the second kind Q„ (A) associated to J are the poly- 
nomials of the first kind for J'l. Matrix Ji also satisfies conditions 
of Theorem [H which yields that there exists a function Fi analytic in 
D\{0} and continuous in D n f/ such that for z G D\{0}, 



(45) Qn[z + - 

z 



zFiiz) 1 



1 



z'^ z^ 



exp 



1-27 



+ 



exp 



1 - 27 



1-27 



27 



as n —7- 00. 
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This and asymptotics ( 133|) of P„ imply that the combination Qn{^) + 
m(A)P„(A) belongs to P for A G C+ and A G C_ only if 

(46) m(^z+^^ = for z eB)\{-l;l). 

It follows that zeros of F in D correspond to eigenvalues of J' outside 
the interval [—2; 2] and hence can only lie on the interval (—1; 1). For 
every z E T (1 U, fHB]) has a limit, 

1 \ Fi(z) 



z J F(z 
This is what we look for, because 



1\ 1 / 1 A F(z)Fi(z) - F(z)FAz) 



zj n \ z J 27ri|F(z)|2 

and the spectrum of J' is purely absolutely continuous on intervals 
(— 2; 2)\{±2 cosw, ±2 cos2ci;} (from the fact that the limit is finite at 
every point of these intervals, P].p^). 
Theorem [1] also gives for z G T D ?7, 



zFAz) „ / ii2{z)n^-^^\ 
— 1 \ 1 — 27 / 

Substituting this expression and the analogous one ( 134|) for P„ into the 
formula for the Wronskian of P and Q (which is constant and equals 
to one), we get after a short calculation: 

= P^(z + -^ Qn+1 Pn+l [Z+ - ]Qn[z+ - 



_ z{F{z)Fi{z) - F{z)Fi{z)) 

1-^2 

where terms o(l) cancel each other, so that the result does not depend 
on n. From this we have: 



F{z)Fi{z)-F{z)Fi{z) = --z, 



z 



which together with fHi|) gives fHSl) . □ 



22 JAN JANAS AND SERGEY SIMONOV 

Corollary 1. Under conditions of Theoreml^ 

n— >oo 

Proof. From ( l35l) . ( l39l) . Lemma [6] and ( HTl) we have: for z G T fl f/, 

(P„+i + i) - zP^ [z + i)) ^"exp (^-i^-'^' 



(P.-,i(. + i)-iP„(.+ i))^exp(^ 



■27 



-27 

F{z) 



+ o(l) as n — 7- oo. 



Therefore 

|P„+i(A) - zPn{\)\ \F{z)\ as n oo, 
and together with f H3|) we obtain the assertion of the corollary. □ 

6. The case 7 g (|; l] 

In this section we formulate the result for the simpler case 7 G (|; l] 
and show how to simplify and modify the proof of the corresponding 
result for 7 G (|; |) to fit this formulation. We will need this as a step 
in proving the asymptotics of the spectral density ([7]). 



Theorem 3. Let {bn}'^=i be given by ([T]) and 

,uj ^nZ and G l^. 



7 € (i; 1 

Then for every z G T\{1, — l,e^*'^, — e^*'^} i/iere exists F{z) such that 
orthonormal polynomials Pn associated to {bn}'^=i have the following 
asymptotics: 



^ I l\ zF{z) 1 zF{z) „ 
Pn(z + -]= ■ — + ■ + o 1 asn^(x. 

z I 1 — z^ z^ z^ — I 



Function F is continuous in T\{1, — l,e^"^, — e^*"^} and does not have 
zeros there. Spectrum of the Jacobi matrix J given by (13]) is purely 
absolutely continuous on (— 2; 2)\{±2 cosw}, and for a. a. A G (—2; 2) 
the spectral density of J' equals: 
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Proof. Consider 7 G l] and return to previous sections. Statement 



of Lemma [T] holds true if we replace U with C\{0, 1, 
the estimate (fT3l) with 



-l,e 



±1^ 



h 



2) 



(^) 



o 



^27 



as n — )■ 00. 



and fi2{z) with zero. In the proof of Lemma [T], we include terms 
of the order O {:^) into the remainder and hence make no use of 
M±4{z),X±4{z),V{z) andrP(^). Condition 

is not needed anymore. System fl20|) becomes 



i 

z 

2 



In statement and proof of Theorem[T]we can also replace U, O + l^n+il) 
andyU2(-2) with correspondingly C\{0, 1, -1,6^"^, -6=^"^}, O + |g„+i|) 
and 0. Most of the estimates that we use there become trivial. The 
same corrections should be applied to the proof of Theorem [2], as well as 
replacing (-2; 2)\{±2 cosw, ±2 cos2a;} with (-2; 2)\{±2 cos w}. □ 
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